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Abstract. We report on the results of an extended optical spectroscopic monitoring campaign on the early-type B supergiant 
HD 64760 (BO. 5 lb) designed to probe the deep-seated origin of spatial wind structure in massive stars. This new study is based 
on high-resolution echelle spectra obtained with the Feros instrument at ESO La Silla. 279 spectra were collected over 10 
nights between February 14 and 24, 2003. From the period analysis of the line-profile variability of the photospheric lines we 
identify three closely spaced periods around 4.810hrs and a splitting of ±3%. The velocity - phase diagrams of the line-profile 
variations for the distinct periods reveal characteristic prograde non-radial pulsation patterns of high order corresponding to 
pulsation modes with / and in in the ran ge 6 - 10. A detai led modeling of the multi-periodic non-radial pulsations with the 
Bruce and Kylie pulsation-model codes iTownsen3l 1 997b() favors either three modes with / = -m and / = 8, 6, 8 or m = -6 
and / = 8, 6, 10 with the second case maintaining the closely spaced periods in the co-rotating frame. The pulsation models 
predict photometric variations of 0.012 - 0.020 mag consistent with the non-detection of any of the spectroscopic periods by 
photometry. The three pulsation modes have periods clearly shorter than the characteristic pulsation time scale and show small 
horizontal velocity fields and hence are identified as p-modes. The beating of the three pulsation modes leads to a retrograde 
beat pattern with two regions of constructive interference diametrically opposite on the stellar surface and a beat period of 
162.8 hrs (6.8 days). This beat pattern is directly observed in the spectroscopic time series of the photospheric lines. The wind- 
sensitive lines display features of enhanced emission, which appear to follow the maxima of the photospheric beat pattern. The 
pulsation models predict for th e two regions normalized flux amplitudes of A = +0.33, -0.28, sufficiently large to raise spiral 
co-rotating interaction regions <Cranmer & Owockll99d) . We further investigate the observed Ha wind-profile variations with 
a simple rotating wind model with wind-density modulations to simulate the effect of possible streak lines originating from the 
localized surface spots created by the NRP beat pattern. It is found that such a simple scenario can explain the time scales and 
some but not all characteristics of the observed Ho- line-profile variations. 
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rotation 



1. Introduction 

The discovery of systematic, patterned variability in the stellar 
winds of luminous hot stars rates as one of the major landmarks 
in recent stellar astrophysics. It is pertinent to our understand- 
ing of wind dynamics and stellar structure. 

Clearly, to make progress understanding these new phe- 
nomena and their impact on the fundamental process of mass- 
loss via stellar winds - and therefore on stellar evolution sce- 
narios - , the mechanism(s) responsible for changing the start- 
ing conditions of the wind and - more specifically - the mecha- 
nisms for dividing the stellar surface into distinct regions must 
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be determined. Furthermore, the underlying physical mecha- 
nisms for coupling the photospheric and wind variations have 
to be worked out in detail. 

The discovery of strictly periodic and sinusoidal modula- 
tions of the Si III, Si iv, and Nv UV resonance lines with pe- 
riods of 1.2 and 2.4 days jPrinia et al.lll995t iFuUerton et alJ 
Il997h and the recent detection of the same periods in the op- 
tical Ha line Jkaufer et al. 2002, KPS 2002 in the following) 
have turned the intrinsically fast rotating early-type B super- 
giant HD 64760 (BO. 5 lb) into a key object in the study of spa- 
tially structured hot star winds and their connection to the stel- 
lar surface: 

KPS 2002 presented the first observational evidence for the 
physical link between the observed photospheric variations and 
the large scale wind structure in HD 64760. Photospheric low- 
order NRP patterns and the constructive or destructive interfer- 
ence of multi-periodic photospheric oscillations were proposed 
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as the most probable source to provide the large-scale pertur- 
bations at the base of the wind needed to build large-scale wind 
structure. The photospheric UV lines of HD 64760 as observed 
intensively with Iue were used to probe for the first time the 
transition zone betw een the deeper photosp here and the base of 
the stellar wind bv iHowarth et alJ ( Il998l) . Unfortunately, the 
study was limited by the comparatively low S /N of the Iue 
spectra. However, this missing link photosphere - wind can 
now be filled observationally with highest quality optical spec- 
troscopic time series, which is the study presented here. 

The optical variabiHty analysis of HD 64760 by KPS 2002 
was based on Heros time series obtained at the ESO50-cm 
telescope and suffered from the comparatively low S /N (~ 100 
at Ha and dramatically dropping towards the blue) and the poor 
time resolution (1 spectrum per night) of the data set. In partic- 
ular, the low quality of the spectra inhibited a detailed study 
of photospheric line-profile variability, which was found to be 
very subtle, i.e., line-profile variations with peak-to-peak am- 
plitudes of less than 0.5% of the continuum level. 

In this work, we present the results from a follow-up obser- 
vation campaign with Feros at the ESO/MPG2.2-m telescope 
at La Silla, which provided us with continuous high-quality op- 
tical spectroscopic time series over more than two rotational 
cycles («i 10 days) with prerequisite high S/N (~ 400) and 
spectral 50000) and time lOmin) resolution. 

The new data set was acquired with the aim to detect and 
to determine the characteristics of the photospheric and wind 
line-profile variations and to further explore the nature of the 
mechanism that couples them together. The scope of this pa- 
per is not to carry out detailed modeling of the observations 
but to provide detailed input for unified hydrodynamic stellar 
wind models, which ultimately need to include pulsation (or 
other physical mechanisms like magnetic fields) as the driver 
for large scale stellar wind structure. 

2. Observations 

The new observations were carried out with the high- 
resolution echelle spectrograph Feros jKaufer et al. 2000) at 
the ESO/MPG 2.2-m telescope at La Silla. A total of 279 spec- 
tra with a resolving power of /? = 48 000 and a wavelength 
coverage from 3600 - 9200 A were collected over 10 nights 
between February 14 and 24, 2003. In a typical night, some 30 
spectra were recorded over time intervals of ~4 hours. Tabled 
provides a more detailed record on the distribution of the indi- 
vidual observations over the different nights. The average ex- 
posure time per spectrum was 400 s but varied between 200 s 
and 820 s to account for the variations of the observing condi- 
tions. A very high average 5/A^-ratio of 390 per spectrum (at 
5400 A) was obtained for this time series. 

Flatfield and wavelength-calibration exposures have been 
obtained with the instrument-internal Halogen and Thorium- 
Argon lamps at the beginning of the respective nights. All spec- 
tra have been reduced semi-automatically with ESO-M idas us- 
ing the dedicated Feros context as described e.g. in lStahl et alJ 
ill999h . All spectra have been reduced to barycentric veloci- 
ties and have been normahzed to the stellar continuum using 



Table 1. Observation log of the Feros observations of 
HD 64760 at the ESO/MPG 2.2-m telescope. 
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the very stable instrument response curve and low-order fits to 
clean stellar continuum points. 

Throughout this paper, all velocities are given with respect 
to the laboratory wavelengths of the corresponding lines of in- 
terest. For this purpose the wavelengths of the respective lines 
have been coiTected by a systemic velocity of v'sy, = -i- 1 8 km s"' 
for HD 64760 as was obtained from the fitting of rotationally 
broadened synthetic line profiles by KPS 2002 (cf. their Fig. 1). 
The uncertainty in Vjys from the fits amounts to +2kms"'. 

Figure^ shows an example of the 35 spectra around the 
Siiii/14553 line as obtained in night4 of the monitoring cam- 
paign. Systematic line-profile variability (LPV) is directly dis- 
cernible and wiU be discussed in detail below. 

3. Spectroscopic time series 

In the following we present the spectroscopic times series as 
so-called dynamic spectra where the time series of spectra is 
displayed in a two-dimensional velocity (wavelength) - time 
image with the intensity levels of the individual spectra repre- 
sented by a grayscale from the low (black) to the high (white) 
intensity cut level. The intensity levels are computed as the dif- 
ference with respect to the average spectrum of the times se- 
ries. Therefore, the discernible absorption (black) and emission 
(white) features in the dynamical spectra are pseudo-absorption 
and pseudo-emission features relative to the time-averaged stel- 
lar spectrum (gray). The average stellar spectrum is constructed 
with the weights of the individual spectra according to their re- 
spective S /A^-ratio. 

In this paper, we present the time series in a collapsed form, 
not showing the full time gap from night to night. The time 
gap remains as a white horizontal bar, while times with no 
data within the dynamical spectra of one night appear as black 
horizontal bars. The seventh night contains only four spectra 
and is completely omitted from the dynamical spectra (but not 
from the analyses below). The dynamical spectra indicate on 
the right axis the observing night and on the left the start and 
end time of the corresponding observations as counted in hours 
from the start time of the exposure of the first spectrum of the 
time series, i.e., t - Ohrs equals to MJD = 52685.0265. In addi- 
tion, all dynamical spectra indicate the projected rotational ve- 
locity of V sin / = +265 km s"' as was determined in KPS 2002 
for HD 64760. 
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Fig. 1. Time series of spectra of the Si in/}4553 line in night 4 of 
our monitoring campaign. The time axis increases towards the 
top and is given in hours since MJD = 52685.0265. The vertical 
bar indicates 5% in intensity with respect to the continuum; 
vsin/ - +265 kms"' is indicated by vertical dashed lines. 



The dynamical spectra of the representative purely photo- 
spheric Si III/14553 line, the photospheric He 1/16678 line, which 
shows some wind emission contribution and the mostly wind- 
sensitive Ha line are shown in Figs. |2 01 and 0] respectively. 
Virtually all available isolated photospheric lines in the spec- 
trum of HD 64760 in the accessible wavelength range from 
3600-9200 A have been examined through their dynamical 
spectra and are found to appear very similar — but very differ- 
ent from Ha. The variations in the photospheric lines are largest 
in the shown Si iii/l4553 and He 1/16678 lines with peak-to-peak 
amplitudes of + 1% relative to the time-averaged spectrum. The 
strong He 1/16678 line does show in addition to the photospheric 
patterns some wind features from Ha, e.g., most prominent, a 
broad pseudo-emission component traveling from the red (day 
3) to the blue side (days 6-9) of the Une. The wind sensitiv- 
ity of the He 1/16678 line is substantiated by a small amount of 
asymmetry in the time-averaged line profile and the existence 
of small emission peaks superimposed to the photospheric line 
profile at -298 km s ' and H-330kms The emission peaks 
reach a maximum intensity of 0.4% above the continuum and 
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Fig. 2. Dynamical velocity - time spectrum of the Sira/i4553 
line. On the right axis the observing night and on the left axis 
the start and end time of the corresponding observations is indi- 
cated in hours since MJD - 52685.0265. The intensities from 
black to white are displayed with + 1 % cut levels. Further, v sin ; 
is indicated by vertical dashed lines. 

are clear mass-loss indicators for a fast rotating wind as was 
demonstrated for Ha in KPS 2002. No such wind features are 
found for the Sim/14553 line supporting its bona fide photo- 
spheric origin. 

The photospheric lines vary strongly within one night, but 
also the variability pattern looks different from night to night 
and changes from periods with low variability (like in nights 
1 and 6) to periods with strong patterns (like in nights 3 and 
9). This observed pattern variability is indicating either a con- 
tribution from non-periodic effects or constructive/destructive 
multi-periodic effects. The pseudo-periods become directly 
discernible if the spectral intensity variations are followed in 
time for a fixed velocity bin (e.g. at zero-velocity) and ap- 
pear to be of the order of the length of one observing night 
(i.e., about 4 to 6 hours). It is important to note here, that 
the length of the observing nights was just sufficient to sam- 
ple the full extent of most of the photospheric variations as 
can e.g. be seen in Figs. and |2l in night 4 where the pat- 
tern at velocity zero evolves from absorption to emission and 
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Fig. 3. Dynamical spectrum of the He 1/16678 line. Displayed 
with + 1 % cut levels, v sin / is indicated by vertical dashed hnes. 

back to absorption within 4.5 hrs. At all times, the photospheric 
lines do clearly show structured variability with alternating pro- 
grade travehng pseudo-absorptions and pseudo-emission fea- 
tures in the dynamical spectra, which are indicative of the 
presence of non-radial pulsations (NRP) in the atmosphere of 
HD 64760. The fact that the slope of the traveling features ap- 
pears to be non-constant and that from night to night periods 
of low and high variability amplitudes can be distinguis hed 
favor the presence of multi-periodic NRPs. iBaadd ( 1 19841) re- 
ported for the first time photospheric line-profile variations in 
highest-quality optical spectra of the Si iii/l4553 and He 1/16678 
lines, which where interpreted as high (\m\ large) and low order 
(\m\ «i 2) NRPs with respective periods of ~ 0. 1 and ~ 0.5 days. 
iHowarth et a l. ( 1998) could detect the 1.2-day period reported 
in the UV wind lines also in the UV photospheric lines of the 
luE Mega data set with a weak signature for a prograde travel- 
ing NRP pattern in the corresponding phase diagrams while 
KPS 2002 report a weak signature of prograde NRPs in the 
optical He 1/14026 line with a 2.4-day period. Only the high 
quality and proper sampling over several contiguous nights as 
achieved with this new spectroscopic data set allows for the first 
time to directly visualize the presence of NRPs in HD 64760. 
The characterization and modeling of the multi-periodic NRPs 
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Fig. 4. Dynamical spectrum of the Ha line. Displayed with 
±3% cut levels, v sin / is indicated by vertical dashed lines. Note 
the extended velocity range of ±800 km s"^ to cover the wind 
features of the line profile. The sharp spectral features are un- 
resolved telluric lines in the Ha region. 

in HD 64760 are the main subject of our analysis in Sects. 0] 
andia 

In contrast to the photospheric lines the dynamical spec- 
trum of the Ha line as shown in Fig.|3]varies very little within 
one night, indicating that predominantly longer time scales of 
the order of days are governing the variability of this line. 
This finding is consistent with the observations in KPS 2002. 
However, so far no evidence could be given for any variabil- 
ity on shorter time scales due to the inappropriate 1-day sam- 
pling rate in the earlier Heros data set used in KPS 2002. The 
new Feros data set provides for the first time high time res- 
olution for the Ha fine together with coverage over several 
nights. The night- to-night variations in Ha show an evolution 
of the pseudo-absorption and emission features. At the same 
time the variations within one night do not reveal the pres- 
ence of any shorter time scales but are consistent with the 
line-profile changes as expected from the night-to-night trends. 
Interestingly an absence of pseudo-absorptions and emissions 
is noticeable in the nights 1 and 5 almost coincident with the 
low amphtudes of variability in the photospheric lines in the 
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Fig. 5. Two-dimensional power spectrum of the Si iii/l4553 line 
around a period of 4.8 hr showing the three significant power 
peaks as identified by the iterative Clean process. The 99.5% 
level of significance of the periodogram is indicated by a (red) 
dashed line. The window function of the data set is shown as 
(blue) dotted line and was shifted to the position and normal- 
ized to the height of the strongest power peak. The strongest 
alias periods can be seen at 4 hr and 6 hr. The three significant 
periods Pi = 4.810hr, P2 = 4.672hr, P3 = 4.967hr are la- 
beled in the right plot, which displays the average power over 
the velocity bin of -v sin / to +v sin /. 



nights 1 and 6. The Ho- line of HD 64760 is strongly affected by 
the mass loss of the supergiant star and the rotation of the wind 
as was demonstrated in KPS 2002 (cf. their Fig. 2). However, 
since their estimate of mass-loss rate of 9 x 10"^ Mo/yr is still 
moderate, it is anticipated to find contributions of the underly- 
ing photospheric profile and its variability on the shorter time 
scales as found for all other photospheric lines (cf. above). This 
is not the case and it is worth noting here, that the absence of 
NRP patterns in the dynamical spectra of Ha will be 'naturally' 
explained by the NRP modeling as carried out in Sect.|5] 



4. Period analysis 

4.1. Photospheric lines 

The qualitative findings from the inspection of the dynamical 
spectra of photospheric and wind lines in Sect.|3will be quanti- 
fied in this section by means of a rigorous time-series analysis. 

As a first step we carried out a two-dimensional period 
analysis on the Feros data set by computing power spectra 
and periodograms for velocity bins of 5 km s"' width across the 
lines of interest. The full frequency range sampled by the data 
set from 1/10 days"' to 1/10 min"' was searched for signifi- 
ca nt periods by application of the Clean algorithm described 
bv ISchwarzenberg-Czernvl ill 996). This algorithm makes use 
of the equivalence of the classical Clean algorithm and l east- 
squar e fits of sine functions to the observed data ( Sch warzl 
|l978) — both performing an iterative deconvolution of the 
window function from the time series. Periodograms and power 
spectra are computed using Lomb-Scargle statistics, which take 
into account for uneven sampling of the data and have a known 
exponential x^C^) probability distribution. Peaks in the peri- 
odograms below a level of significance of 99.5% probability 
were rejected from the analysis. For a more detailed de scrip- 
tion of this Clean implementation cf. lKaufer et alJ lll997l) . 

Figure |5] shows the section of the two-dimensional power 
spectrum where the most significant periods are found in the 
photospheric lines, i.e., close to a period of 4.8 hr. Application 
of the iterative Clean algorithm as described above confirms 
that three significant periods with Pi =4.810hr,P2 = 4.672 hr, 
P3 - 4.967 hr are present. The three periods are sorted by 
their associated power and significance with Pi and P2 be- 
ing of comparable level but containing about twice as much 
power than P3. The power within each detected period is not 
evenly distributed over the velocities but strongly peaked at 
zero velocity indicating stronger variability at the line center 
compared the wings of the line. Since the actually measured 
periods show some variations of the order of 1.25% as function 
of the velocity within the line profile, the periods given above 
have been measured from the position of the power peaks aver- 
aged over the full velocity interval -v sin ; to +v sin i.e., -265 
to -H265kms"' (cf. Fig.|51 right plot). 

To assure ourselves of the detected multi-periodicity of the 
photospheric variations we have to exclude that they are pos- 
sibly caused by the complex window function of our observa- 
tions. This is implicitly guaranteed by the Clean method, which 
we have used to identify the periods. However, we have veri- 
fied the multi-periodicity in addition as part of the pulsation 
modeling of the variability, which will be described in Sect.|5] 

The period analysis confirms the multi-periodicity of the 
photospheric line-profile variations in HD 64760 as had already 
been suspected from Fig. |2l It should be noted that the three 
periods are closely spaced in time with +3% separation and 
two periods show similar power, which are the prerequisites 
for the possible existence of beat periods on time scales much 
longer than the underlying periods themselves. 

The two-dimensional velocity - phase diagrams of the pho- 
tospheric Si III/14553 line constructed with the three identified 
periods are shown in Fig.|6land clearly show the characteristic 
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Fig. 6. Dynamical phase spectra of the Si iii/14553 hne. The observed spectra are phased with Pi - 4.810hr (left), P2 - 4.672 hr 
(center), Pj, - 4.967 hr (right). 25 phase bins were used for the phase interval 0.0 - 1.0. All spectra are displayed with +0.3% 
cut levels. The measured typical acceleration of the features over the line center of {dv/d(f)) = - 175 km s" '/cycle is indicated by 
a dot-dashed lines. 



traveling features of three distinct prograde NRP modes. While 
the associated pulsation modes of Pi and P3 appear to be very 
similar, the pulsation mode of P2 is clearly different as can be 
seem from the different slope of the NRP pattern. 

Alternatively to the identification of the origin of the line- 
profile variations (LPV) as non-radial pulsation patterns, sim- 
ilar LPVs could be produced by the rotational modulation 
by regular surface structures. For Pi a crossing time of an 
individual surface feature from -vsin; to H-vsin; is deter- 
mined to be 13hrs (according to a measured acceleration of 
200km s"' /4.81 hrs at the line center). This would corresponds 
to a stellar rotation period of Prot = 26 hrs. The breakup rota- 
tion velocity for HD 64760, however, is estimated to Prot.break = 
64 hrs assuming a stellar radius of P = 22 R0 and mass of 
M = 2OM0 (cf. KPS 2002). From this simple estimate it must 
be concluded that the observed LPVs cannot be explained 
by co-rotating surface structures. Further, the observation of 
closely spaced multiple periods is not compatible with rota- 
tional modulation where only the stellar rotation period or in- 
teger fractions of the rotation period can be observed. 

Therefore, we have to conclude that the most likely source 
for the observed LPVs are multi -periodic non-radial pulsations. 
In Sect. |5] we will further substantiate this hypothesis by de- 
tailed modeling of the NRPs. 

4.2. The wind-sensitive Ha line 

The Ha line has been examined for periodicities in the same 
way as the photospheric lines using the iterative Clean algo- 
rithm described above. The power spectrum of Ha is very com- 
plex in particular in the frequency regime of < 1 day ' where it 
is hardly possible to identify the most significant frequencies to 
which the Clean procedure shall be applied. Therefore, we had 
to closely inspect selected frequency ranges for which period- 
icities could be expected. First, we tried to verify the presence 



of the three photospheric pulsation periods close to 4.8 hrs. But 
no significant power could be found as expected from the lack 
of fast moving features in the dynamical spectrum of Ha. 

The short extent in time of the new Feros dataset(10 nights 
with several spectra within a few hours of each night) is not 
very well-suited to search for the continued presence of the 2.4- 
day wind modulation period as it was found in the Iue Mega 
data set (15 days with continuous sampling) and confirmed 
in Ha with the Heros data set (100 nights with 1 spectrum 
per night). Nevertheless, the CLEANed two-dimensional power 
spectra show significant power at a period of 2.4 days. If the 
Ha data sets of the Heros (1996) and the Feros (2003) are both 
phased with the established period close to 2.4 days, a strik- 
ingly similar variability pattern is found for both optical data 
sets as can be seen in Fig. Both data sets have been treated 
in exactly the same way by sorting the spectra in 25 phase bins 
(for phase 0.0 - 1.0) according to a period of 2.39097 days. 
Phase 0.0 is defined as the time of the observation of the first 
spectrum of the 2003 time series, i.e., MJD = 52685.0265. The 
apparent accuracy of the period comes from the requirement to 
match the 1996 and the 2003 phase diagrams. However, this re- 
quirement is fulfilled by a set of discrete periods. Using phase- 
dispersion minimization of the combined 1996+2003 data set, 
we have identified 20 such possible periods in the range from 
2.37-2.41 days. Both time series have been resampled to a ve- 
locity resolution of 5 km s" ' and are displayed with ± 1 % inten- 
sity cut levels. The respectively subtracted time-averaged Ha 
spectra from the 1996 data set (cf. KPS 2002, their Fig. 1) and 
the 2003 data set to compute the shown residuals are indis- 
tinguishable within the limitations of available S /A^-ratio and 
achievable accuracy in the normalization of the stellar contin- 
uum. The variability pattern of the 2.4-day phased Ha line in 
Fig.0is very characteristic with 'X' or '8'-shaped pattern for 
the pseudo-emission features with the 'crossing point' close to 
phase 0.20 (or night 6 in the 2003 data set). KPS 2002 described 
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Fig. 7. Dynamical spectra of the Ha line phased with a period of 2.39097 days (cf. text). 25 phase bins were used for the phase 
interval 0.0 - 1.0. Left 1996 Heros data set, right 2003 Feros data set. Both dynamical spectra are displayed with +1% cut levels. 



the 2.4-day Ha variability in detail and tried to disentangle the 
different contributing components at distinct velocities in the 
line profile. All the features highlighted there can also be identi- 
fied in the new Feros data set — however, the features appear to 
be more coherently traveling across the line profile. The resem- 
blance of the two dynamical Ha-spectra in Fig.0indicates that 
2.4-day period has now been stable for at least eight years, i.e., 
the time span between the 1995 Iue, the 1996 Herds and the 
2003 Feros campaigns. Such a long-term presence of a single 
modulation period requires a stable clock in the system. Stellar 
rot ation has been favored by aU previous works on HD 64760, 
cf. iFullerton e t aT ('1997") and KPS 2002. In this work we will 
further investigate the possibility that stellar oscillations could 
act as stable clock in HD 64760, too. This possibility will be 
discussed in context with our pulsation and wind modeling in 
Sects. Eland 



Table 2. Bruce/Kylie model parameters 



Model Parameters 
Tpoic 29 000 K temperature at pole 

M 20 Mg total mass 

iJpoie 18Ro radius at pole 

Vequator 265 km s"' rotation velocity at equator 

( 90 deg inclination 

Derived Parameters 
/^equator 21.6 R© radlus at cquator 

log^poie 3.23 gravity at pole 

log^fequatoi 2.85 gravity at cquator 

requator 23 300 K temperature at equator 

Observable Parameters 
T^i 98.9 hr (4.1 2d) rotation period 

L,ot 1.55xlO^Lo luminosity 

Teff 24 600 K effective temperature 



5. NRP line-profile modeling 

The period analysis of the photospheric lines has revealed the 
presence of multi-periodic variabiUty with three dominant pe- 
riods. We have identified non-radial pulsations (NRPs) as the 
most likely origin of this variability and proceed now to con- 
firm this hypothesis with detailed modeling of the LPVs. 

In the following we first concentrate on the determination 
of the pulsation quantum numbers / and m of the NRP modes 
associated with the three pulsation periods. Only detailed mod- 
eling will allow us to test the consistency of the very complex 
observed LPVs with multi-periodic NRPs. However, because 



of this high complexity of the observations on one hand side 
and the very large parameter space on the other, most likely, an 
exact identification of the pulsation parameters is not possible 
with our observational data set. 

We u se the NRP modelin g codes Bruce and Kylie (Version 
2.84) by iTownsenj ( Il997bl) to examine the LPVs observed 
in the photospheric lines of HD 64760. Bruce calculates the 
pulsational perturbations to the rotating stellar surface and in- 
cludes pulsational velocity fields in all three spatial directions 
(r, (p, 6), temperature fluctuations, variations in the visible sur- 
face area and the variations in the viewing angle of a given 
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surface element. Kylie derives in a second step the observable 
quantities from the perturbation fields provided by Bruce by 
computing and integrating the observables like the observer- 
directed spectrum or flux of each surface element individually 
from an input grid of stationary atmosphere models. 

The intrinsic line profiles and spectral energy distributions, 
which are required as input to Kylie for the production of 
realistic non-radial pulsation LPVs and light curves, respec- 
tively, were computed on the basis of a n extended grid of 
Atlas 9 model atmospheres ( Kurucz'1993') and subsequent Bht 
spectrum synthesis (Baschek et al. 1966) for solar abundances 
and a microturbulence of ^mkro = 2kms"'. The grid extends 
from 15 000 to 31 000 K in steps of lOOOK in effective tem- 
perature Tgff and from 2.6 to 3.5 in steps of 0.1 in gravity 
logg. For the hotter atmospheres (25 000 K and higher), the 
Atlas 9 models with the lowest log g values do not converge 
anymore. Fortunately, these grid points are not required in the 
Kylie line-profile computations, since for a fast rotating star 
like HD 64760 the hottest regions are located at the stellar pole 
where the log g is high because of the rotational flattening of the 
star. The lower temperatures occur close to the equator, where 
also log g is found to be lower. 

The stellar input parameters for Bruce were chosen in such 
a way to reproduce the observable parameters of HD 64760 as 
given in KPS2002 (their Tab. 1): a polar radius of 18 Rq, a 
polar temperature of 29 000 K, a stellar mass of 20 M ©, an 
equatorial rotation velocity of 265 km sec"' (65% of the critical 
rotational velocity of 415 km sec"') and an inclination of 90°. 
These input values result in an equatorial radius of 21.6 Rq, 
logg values of 3.23 and 2.85 at the pole and equator, respec- 
tively, an equatorial temperature of 23 300 K, and a rotation 
period of 98.9 hours (4. 12 days). The total luminosity and effec- 
tive temperature, as seen from the observer, are 1.55 x 10^ 
and 24 600 K. The model parameters related to the stationary 
stellar model are summarized in Tab.|2 

5.1. Pulsation mode parameters 

For the modeling of n pulsation modes, Bruce needs the pulsa- 
tion periods P„ as measured in the observer's system, the pul- 
sation amplitudes A„, the relative pulsation phases (p„ and the 
pulsation quantum numbers m„ and l„. 

The pulsation quantum number m is equal to the number of 
meridian circles on the stellar surface and can be estimated via 
the acceleration of the traveling features across the line center 
From the phase diagrams on Fig. |6lwe measure an accelera- 
tion of {dvld(p) - -175 + 50kms"'/cycle. This translates with 
the purely geometrical relation m„(dv/d<p)n/2 = 2|vsin/| into a 
possible range of values for m of -8.5 < m < -4.7. 

While the m-value can be estimated with a reasonable ac- 
curacy by this method, the Z-value remains ill-defined and 
can only be determined by detailed modeling of the LPVs. 
However, the presence of multiple pulsation periods and the 
high but still finite S /A^-ratio of the data will naturally limit the 
achievable accuracy of the /-value determination. 

The amplitudes and relative phases of the pulsation modes 
have been chosen so that the best possible match to the obser- 



Table 3. Bruce/Kylie model parameters 



Pulsation Modes 



n 


P„ [hr] 


f*iixorot [hr] 


<l>„ [deg] 


/„ 


m„ 


A„ [kms '] 


Model 1 


1 


4.810 


7.872 


180 


8 


-8 


5 


2 


4.672 


6.520 


216 


6 


-6 


5 


3 


4.967 


8.302 


97 


8 


-8 


3 


Model 2 


1 


4.810 


6.791 


180 


8 


-6 


10 


2 


4.672 


6.520 


216 


6 


-6 


5 


3 


4.967 


7.109 


97 


10 


-6 


10 


Model 3 


1 


4.810 


7.872 


180 


6 


-6 


5 


2 


4.672 


6.520 


216 


8 


-8 


5 


3 


4.967 


8.302 


97 


6 


-6 


3 



vations is achieved: pulsation amplitudes of 3 - 10km s"' are 
needed to produce LPVs of the observed amplitudes of +1% in 
the Si III/14553 line; the relative phases have been adjusted to 
match the times of strongest variations in the observed time se- 
ries, i.e., the times of constructive interference of the pulsation 
modes. 

Our adopted pulsation mode parameters are summarized in 
Table |5] Two fundamental types of possible pulsation modes 
have been investigated, / - -m (Models 1 and 3) and / -m 
(Model 2) and are discussed below. 

The measured and derived stellar and pulsation-mode pa- 
rameters in Tabs. |2l and |3] allow us to identify the expected 
pulsation-mode characteristics. The dimensionless pulsation 
frequency in the co-rotating frame Wcorot is defined as 



(27r/Peorot)'-R? 



pole 



GM 



(1) 



and reaches values of Wcorot ~ 30 for the observed pulsation 
parameters of HD 64760. This places the three observed pul- 
sation modes in the regim e of the pressure- or p-mode pul- 
sations jTownseridll 1 997ah . This finding is consistent with the 
observed concentration of power at the line center in the two- 
dimensional power spectrum (Fig.|5} indicating a strong prefer- 
ence of radial over horizontal velocity fields, which is the most 
direct characteristics of p-mode pulsations. In the NRP models, 
the ratio of the horizontal to the radial velocity amplitudes k is 
described as A: = l/w^^j.^, and will provide a good representa- 
tion of the power distribution over the line profiles (cf. Fig. 1101 
below). 



5.2. The case I = 



-m 



The analysis of the phase diagrams above provides a reason- 
able estimate of the m-value of each pulsation period but the 
degree / remains ill-defined, m can assume values from -I to 
+1, i.e., the lowest acceptable value of / in any case is -m. For a 
given choice of /, there will be I - \m\ nodes along any longitu- 
dinal line, and consequently, the stellar surface will be divided 
into I - \m\ + 1 sections in latitude. With an increasing number 
of zones on the stellar surface the mutual cancellation of the 
pulsation effects of neighboring zones will increase. Therefore, 
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Fig. 8. Computed time series for the Siiii/14553 line accord- 
ing to NRP Model 1. Note the strong changes from night to 
night, which are due to the superposition of modes with dif- 
ferent m values. The choice of the m- values for the different 
periods leads in this model to a retrograde propagation of the 
beat pattern - consistent with the observations in Figs.|2land|3] 
Displayed with 1 % cut levels. 

to produce the maximum observable effects, the smallest num- 
ber of surface zones should be chosen, i.e., the case I = -m. 
This is the reason, why in many works on NRPs, this choice 
is preferred. However, usually the same amplitudes in the ob- 
servables but with different characteristics can be achieved with 
higher amplitudes of the modulations in velocity or tempera- 
ture in the case / 9^ m. 

We assume in our first test model (Model 1) / = -m and 
; = 6 for the P2 = 4.672 hr and / = 8 for the Pj = 4.810hr 
P3 = 4.968 hr periods. For pulsations with different m-values, 
an interesting pattern for the resulting total pattern is produced: 
several temperature and velocity maxima emerge along the 
equator from the constructive or destructive interference of the 
pulsation modes. Depending on the period difference of the m- 
modes, this beat pattern moves across the stellar disk in the 
direction of the rotation (cf. our Model 3, Fig. I12> . but also 
against this direction (this Model 1, Fig.|8j. The relevant pe- 
riod is the beat period of the two modes. 



— Lclml +ijGlriL 
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Fig. 10. Two-dimensional power spectrum of the Sim/i4553 
line around a period of 4.8 hr for Model 1. The power distri- 
bution in velocity and frequency of the observations (Fig. |5} 
is reproduced very well by this model. The concentration of 
power at the line center is due to the p-mode characteristics of 
the pulsation modes. 

For a best possible comparison with the observations, we 
computed the NRP models with exactly the same sampling pat- 
tern as during the observations. Dynamical spectra of the com- 
puted time series of the Sini/14553 line and the photospheric 
part of the Ho- line for Model 1 are shown in Figs. Islandll II 
respectively. The model for the Sim/i4553 line reproduces the 
strong changes from night to night in the observed variabil- 
ity pattern. While the dominating non-radial pulsation features 
travel prograde over the line profile within a few hours, the con- 
structive interference of the Pi and P2 pulsation modes travels 
retrograde over the line profile with a period determined by the 
beat period of the two dominant pulsation modes: 

1 _ 1 1 _ 1 1 _ 1 

~ P^'Vi ~ 4.672 hr ~ 4.810hr ~ 162.8 hr ^'^^ 

The maximum amplitude of the constructive interference is in- 
deed observed on the nights 3 and 9 corresponding to the beat 
period of 162.8 hr or 6.8 days. Further, the observed phases 
with low variability in nights 1 and 6 (cf. Figs.|2]|3|l are quite 
well reproduced by this model with three beating periods. It is 
important to note here the directly observable secondary effects 
on the stellar wind as seen in the superimposed wind emission 
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Fig. 9. Dynamical phase spectra for the the Siiii/14553 line as computed for Model 1. The spectra are phased with Pi = 4.810hr 
(1 = S,m = -8), P2 = 4.672 hr (I = 6,m = -6), P3 = 4.967 hr (/ = 8,m = -8). 25 phase bins were used for the phase range 
0.0 - 1.0. All spectra are displayed with ±0.3% cut levels. 



■Triours -vajfji':. ^ K'Gir: night TiOurs ■'v-siri +'-sir: 1 ■■^lil. 




Fig. 11. Computed time series for the photospheric contribution Fig. 12. Computed time series for the Si iii/l4553 line according 

of the Ha line according to Model 1 . The expected variations in to Model 3. The fact that the larger m-values has been assigned 

the line due to NRPs are very small consistent with the obser- to the longer pulsation period leads to a prograde travehng beat 

vations where the pulsation periods could not be detected. The pattern, which is not consistent to the observations. Displayed 

strong variations at the right side are due to the neighboring with 1 % cut levels. 
Cia/l6578,6582Unes. 
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in the|3line and - even more pronounced - in the Ha line. The 
emission feature travels across the line profile from red to blue 
and seems to follow the retrograde pattern and is strongest in 
the nights 2 to 4 and 6 to 9 . 

Figures |5] and show the phase diagrams and the two- 
dimensional power spectrum as computed from Model 1 and 
have to be compared to the observations in Figs. |6l and |5j re- 
spectively. An excellent agreement with the observations is 
achieved with this multi-periodic NRP model, including the 
p-mode characteristics of the pulsation modes. On the other 
hand, Model 3, which produces a prograde traveling beat pat- 
tern (cf. Fig.ll2> does not represent the observations very well 
and must be excluded. 

The weakness of the predicted variations in the Ha line as 
shown in Fig. ^2 (displayed with the same intensity cut lev- 
els as the Si iii/14553 line in Fig. |S} is due to the intrinsically 
(Stark-) broadened photospheric Ha-line profile, which tends 
to smear out the NRP-induced line-profile variations over the 
extended line profile. This at first glance surprising result of 
the NRP modeling is consistent with the non-detection of any 
of the three NRP pulsation periods in the observed Ha time 
series. 

On the other hand, the finding of NRP beat periods in the 
order of several days is of great importance for the investigation 
of the possible connections of (fast) photospheric variability on 
time scales of hours to (slow) wind variability on time scales of 
days. The possible implications will be discussed extensively 
below in Sect.|6lwhere we will try to understand and model the 
effects of slowly varying photospheric surface structures on the 
stellar wind. 

5.3. The case I i= -m 

The observed periods 2,^*3 are closely spaced in time in 
the observer's reference frame. However, since the transforma- 
tion of the periods to the co-rotating reference frame 



'rot 

only depends on the pulsation quantum number m, the pulsa- 
tion periods are actually very different in the co-rotating frame 
for the case / = -m as discussed above in Model 1. The ob- 
served close spacing of the observed periods would be in this 
case merely an unlikely coincidence. 

Therefore, we tested with Model 2 the case of / -m but 
with m identical for all three pulsation modes since in this case 
the close spacing in the observer's reference frame is naturally 
explained by a close spacing in the co-rotating reference frame. 

In order to maintain the interesting beat patterns, which ap- 
pear to be essential to explain the observed NRPs in the pho- 
tospheric lines, we have to assume different /-values for the 
individual modes. In Fig. we show the resulting line pro- 
file variations for one such case, with Z-values of 10, 8 and 6 
and an m-value of -6 (Model 2 in Tab. O. It can be seen that 
qualitatively similar changes can be produced as for the case 
of I = -m and different /, m-values. With slightly larger veloc- 
ity amplitudes of the pulsations the same amplitude of profile 
variations and the same power distribution over the hne profile 



riDOBl 2 
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Velocity [Km/sl 

Fig. 13. Computed time series for the line Si iii/l4553 according 
to Model 2. Note the strong changes from night to night, which 
are due to the superposition of modes with different / values. 
Displayed with 1% cut levels. 



is produced. The increased velocity amplitudes are needed to 
compensate for the cancellation effects from surface zones of 
opposite phase. However, despite the fact that the two inves- 
tigated cases with properly chosen model parameters lead to 
practically indistinguishable phase and power spectra, the ac- 
tual resulting surface patterns are clearly more complex in the 
I + -m case and in particular display strong differences be- 
tween pole and equator within the same meridional sector of 
the star 

5.4. Multi-line analysis 

One possibility to discriminate between I - \m\ and / 9^ \m\ pul- 
sation modes is a multi-line analysis of the NRPs. The strong 
variations of temperature and gravity from pole to equator in 
a fast rotating and pulsating star like HD 64760 provide addi- 
tional diagnostics by the analysis of spectral lines with different 
temperature and gravity sensitivity. With polar/equator temper- 
atures of 29 000/23 300 K, the temperature sensitive He ii/14686 
hne is of particular interest in this case. 
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Fig. 14. Dynamical spectrum of the He ii, O ii spectral region 
computed according to the NRP Model 2 assuming adiabatic 
temperature effects (cf. observed data in Fig. I15> . +v sin / for 
the O II (left) and the He ii (right) lines are indicated. 

FigureOshows the synthetic profiles for the the Hen line 
and the neighboring O ii/l4676 line according to the parameters 
of our I + \m\ NRP Model 2. While the global appearance of 
the two lines is the same, distinct differences can be seen in the 
details of the NRP patterns, in particular in the smaller slope 
and the lesser extent in velocities of the features (i.e., they do 
not reach +v sin /). These characteristics are as expected for the 
Hen line, which is predominantly formed closer to the stel- 
lar pole at higher temperatures: (i) the number of features at 
a given time remains the same, (ii) v sin / becomes apparently 
smaller for that line since at the polar caps only smaller pro- 
jected rotation velocities are seen, (iii) the slope of the travel- 
ing features has to become smaller since the pulsation periods 
have not changed. 

The observed dynamical spectrum as shown in Fig.llSlcor- 
responds well for the O ii line as expected due to its similarities 
to the Si III line for which the model was constructed before and 
a good match with the observations was achieved. However, the 
model does not match for the He n line. In particular the ampli- 
tudes of the variations of the He ii line are observed to be very 
small and hardly discernible with the given S/N while in the 
model the amplitudes are predicted to be even larger than for 
the On line. 

The fact that the models predict too large variations for the 
temperature sensitive Hen line while less temperature sensi- 
tive lines are well represented indicates that the Bruce mod- 
els overestimate the pulsation-induced temperature variations. 



Fig. 15. Dynamical spectrum of the Feros observations of the 
Hen, On spectral region. This dynamical spectrum has to be 
compared to the model computations as shown in Figs. 1141 and 
[T6l +v sin / for the O n (left) and the He n (right) lines are indi- 
cated. 

The models assume adiabatic temperature variations accord- 
ing to .Buta_& S mith (1979). A good agreement between ob- 
servation and model can be achieved if the temperature effects 
are reduced to 50% of the adiabatic value, which leads to a 
reduction of the amplitude of the variations by a factor of 5. 
The corresp onding dynamical spe ctrum is shown in Fig. [T6l 
Interestinglv. lButa & Sm ith* ('1979') found already in their orig- 
inal work on the non-radial pulsator 53 Persei that the assump- 
tion of adiabaticity led to an overestimation of the pulsation- 
induced photometric variations by one magnitude. From our 
observations we have to conclude that the assumption of adi- 
abaticity also breaks down for realistic pulsation models of 
HD 64760. 

5.5. Photometry 

Another possibility to discriminate between I - \m\ and / \m\ 
pulsation modes is to compare the observed photometric vari- 
ations with the photometric variations as predicted from NRP 
models. 

Figure shows the results from our Bruce/Kylie NRP 
simulations for the Johnson V filter for the two studied cases 
of NRP modes (i.e.. Model 1 and 2 as defined in Tab. O. For 
both cases, the maximum (peak-to-peak) amplitude of the vari- 
ations is 20mmag which is consistent with the photometric 
variations of up to 30 mmag as frequently found in OB super- 



A. Kaufer et al.; Multi-periodic photosplieric pulsations and connected wind structures in HD 64760 



13 



□ 1I\4S7S HEn?L4^S 




JSETO ^sb:j 

l^ue length [1] 



Fig. 16. Dynamical spectrum of the He ii, O ii spectral region 
as computed from Model 2 but with the pulsation-induced tem- 
perature effects reduced to 50% of the adiabatic value, (cf. ob- 
served data in Fig. I15>. +vsini for the On (left) and the Hen 
(right) Unes are indicated. 



giants the temperature effects in the NRP 

models are reduced to 50% of the adiabatic value as suggested 
by our analysis in Sect. 15.41 this leads to a reduction of the max- 
imum amplitude of the photometric variations to 12 mmag. 

While the high-frequency variations in the light curves are 
directly related to the dominant individual NRP modes with 
periods close to 4.8 hours, the overall evolution and the shape 
of the envelope of the light curves is governed by the beat- 
ing of the pulsation modes. In particular the clearly differ- 
ent characteristics of the light curves at times of constructive 
(f 60, 220hrs) and destructive interference (f « 160 hrs) 
is of high potential to discriminate between the two studied 
cases: while in the case of / = \m\ (Model 1) the amplitude 
of variations is minimal at times of constructive interference, 
the amplitude is maximal in the case of / \m\ (Model 2) 
and vice versa at times of destructive interference. Further, 
while Model 1 predicts an increase of the average brightness 
by 10 mmag at times of constructive interference, the variations 
of the average brightness according to Model 2 are reduced to 
some 5 mmag without clear relation to the evolution of the beat 
pattern. 

1 1 1 [ [ \ i< V 1 1 photometry han Leeuw en et all Il997 i) of 
HD 64760 (151 measurements over 550 days) displays disper- 
sions of 7, 26, and 22 mmag in the Hp, Bj and Vj passbands, 
respectively. The measured variations are comparable with our 
model predictions. However, no significant periodic signals 
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Fig. 17. Photometric variations in the Johnson V band as pre- 
dicted by NRP Model 1 (upper light curve) and NRP Model 2 
(bottom light curve) for the full 10-day observing run and 
with a time sampling of 30 minutes. The magnitude scale of 
the model has been adjusted to match t he V magnitude of 
HD 64760 from the literature (V = 4.24, iHoffleit & Jaschekl 
The synthetic lightcurve for Model 1 has been shifted 
by -0.025 mag for better display. 



can be found in the data sets. Unfortunately - to the authors' 
knowledge - no other suited photometric data sets are available 
for HD 64760 to verify the predictions of the NRP models. 
The required photometric accuracy of a few mmags appears 
technically feasible for a dedicated photometric monitoring 
campaign, which could put additional rigid constraints on the 
NRP modeling, in particular on the selection of the pulsation 
mode parameters. 

6. Stellar-wind modeling 

One of the key questions for understanding the nature and ori- 
gin of the wind variability in hot stars is its connection to the 
stellar photosphere. It is widely assumed, that the winds are 
modulated by a mechanism related to the photospheric rotation, 
presumably patches on the stellar surfaces produced either by 
non-radial pulsation (NRP) patterns or magnetic surface struc- 
tures ("spots"). The spatial structures in the photospheres lo- 
cally change the lower boundary condition of the stellar wind, 
which in turn causes localized structural changes in the wind. 
These structural changes modulate the observed stellar wind 
profiles as they are dragged through the line of sight to a dis- 
tant observer by underlying rotation of the star 

In KPS 2002 we have introduced a simple 1 .5D-model 
(dubbed Rotbalmer) for a differentially rotating stellar wind to 
model the double emission Ha profile of HD 64760. The model 
is based on the work and description of Petrenz & Puis ( 1996) 
and allows a correct handling of the rotationally twisted reso- 
nance zones in the line formation of wind lines like Y\a. The 
model makes the simplified assumption of a spherically sym- 
metric density stratification throughout the wind. 

To estimate the eff'ects of a geometrically structured rotat- 
ing wind, we have introduced in the wind model density en- 
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Table 4. Rotbalmer model parameters 



1( 



Model Parameters 



T 

J eff 




effective temperature 


p 

^^-equator 


1 A P 


radius at equator 


^equator 


9fiS km t;"' 


VCHJL.liy dL CUU-d-lLfi 


I 


90deg 


inclination 




1500 km S-' 


wind terminal velocity 


^min 


lOkms-' 


wind start velocity 


P 


0.8 


wind acceleration parameter 


M 


9 X 10-^ Mo/yr 


mass-loss rate 


^ He 


0.10 


helium abundance n(He)/n(H) 


iHe 


2 


free electrons per Helium atom 




Streak-line Parameters 






Model 1 


«s 


2 


Number of streak lines 


Pc 


1.3 


Density contrast w.r.t. ambient 


0s 


50deg 


Angular extension (2 arctan 0.5) 


Vs 


265kms-' 


Foot-point rotation velocity 


Model 2 


Is 


2 


Number of streak lines 


Pc 


1.3 


Density contrast w.r.t. ambient 


0s 


50deg 


Angular extension (2 arctan 0.5) 


Vs 


-80.5 km S-' 


Foot-point rotation velocity 



/ 




Fig. 18. Streak-line geometry for Viot - H-265kms ' and 
v's = -)-265kms"' (Model 1) and v^ot = -1-265 km s"' and 
v's = -80.5 km s"' (Model 2). The observer looks at the star 
from the right at model phase (p - 0. 



hancements with respect to the local am bient wind fo l lowin g 
a streak-line geometry as described by iFullerton et alJ (Il997l) . 
The streak line is characterized by the density contrast pc with 
respect to the ambient wind, the angular extension 0, on the 
stellar surface and the rotation velocity Vj of the foot point. 
Usually, Vs is set equal to the rotation velocity Vjot of the star 
corresponding to a streak line emerging from the stellar surface 
and therefore locked to the stellar rotation. 



6. 1 . The case Vs = v^t 

The most natural clock for a strictly periodic wind modula- 
tion is provided by the stellar rotation. For HD 64760 the upper 
limit for the stellar rotation period of 4.2 days is within the ex- 
pected errors close to twice the 2.4-day modulation period. If 
this estimate is correct, the modulating structure on the stel- 
lar surface must be also persistent over such a time scale and 
geometrically divide the stellar surface in two structures with 
180 degree azimuthal separation. 

Our first Rotbalmer wind test model (Model 1) tries to 
simulate this situation with the presence of Ws =2 diametri- 
cally opposed streak lines in the equatorial plane and with their 
foot points rotating at the same speed as the stellar surface, i.e., 
Vs - Vrot : = Vequatoi sin /. The complete set of Rotbalmer model 
parameters used in the model is given in Tab.^and is basically 
identical to the parameter set used in KPS 2002 except for the 
additional parameters to describe the streak lines. The geome- 
try of the emerging streak lines is shown in Fig. 1 181 

Figure shows the dynamical Ho- spectrum as produced 
by this model with two co-rotating streak lines. We have again 
computed the model spectra with the same timings as the ob- 
servations to ease their comparison. If the model is compared 
to the observations in Fig.|4]it is found that neither the veloc- 



hours 
225,2 



"i ight 



f- 1 1 1 1 1 ' ' 1 ' - 


w \ ■ 

It I 
f 'a 


'J 

i; 

\— 





10 




Fig. 19. Ha Rotbalmer Model 1 with Vj 
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Fig. 20. Ha Rotbalmer Model 2 with Vs = -80.5 km s ' corre- 
sponding to a beat period of 162.8 hrs. 



ity evolution of the pseudo-emission and absorption features 
(the 'X' or '8' shaped pattern) is well represented nor the time 
scales of the variations. The dominant time scale is due to the 
presence of two streak lines equal to half the stellar rotation 
period, i.e., 4.2/2 = 2.1 days - the model had been designed 
to exactly match this condition. The model further predicts 
Ha line-profile variations with blue- and redwards accelerat- 
ing features as expected for two symmetric streak lines rotating 
through the line of sight of the observer The features are seen 
from zero velocity up to the terminal velocities of the wind. 
This result is c omparable to the ^ Puppis (04 1(n)f) models by 
lHarriesl ( l2000l) who investigated the effects on the Ha line pro- 
file of a narrow and high-contrast one-armed spiral in the ro- 
tating wind of thi s prototypical O-type star Indeed, we have 
used the model bv lHarriesI ll2QO0l) to verify the functionality of 
our simpler Rotbalmer code and can report that we find a good 
agreement in the modeled LPVs. 



6.2. The case Vs i= Vmt 

Our second Ha rotating structured wind model assumes again 
streak lines spaced at equal angles along the equator However, 
in contrast to the previous case, we do not assume that the foot 



points of the streak lines are fixed on the star and follow its rota- 
tion. In contrast, we assume that the streak lines follow the pul- 
sation pattern, as resulting from the superposition of the three 
detected NRP modes. This pattern does not follow the rotation, 
but has a pattern speed slower than the rotation or - as in our 
case - can even move (in the observer's reference frame) in the 
opposite direction as the rotation. While this does not change 
the stream line of an individual particle, the streak lines can 
change significantly (cf. below). The number of streak lines 
depends on the pulsation pattern. In the case of I = -m, two 
interfering pattern with different m produce nii - OT2 maxima 
along the equator If I + -m and m equal for both modes, m 
maxima along the equator appear For arbitrary values of / and 
m (or more than two interfering modes), the interference pat- 
tern can become very complex. Therefore, we investigate here 
the case / - -m (corresponding to our NRP Model 1), which is 
more attractive as lower boundary condition for a wind model 
with equatorial streak lines because it favors larger pulsation 
(beat-) amplitudes at the stellar equator This second Ha rotat- 
ing structured wind model assumes again - 2 streak lines 
diametrically opposite at the equator corresponding to the 2 
maxima in the beat amplitudes as created for mi - 1112 - -2 

The period in which the beat pattern moves across the star 
is related to the beat period. Therefore, periods significantly 
longer than the individual pulsation periods can result. The beat 
period is not related to the rotation period of the star, but only to 
the pulsation periods. If we consider the two strongest pulsation 
modes Pi and P2 only, we obtain Pbeat = 162.8 hr (cf. Eqn.|2}. 
The apparent rotation velocity of the maxima of the beat pattern 
is then computed as: 



Vi-ot('«l - '«2) 



_Pro 

Pbe; 



265 km s"' (-8 - -6) 
-80.5 km s"' 



98.9 hr 
162.8 hr 



The geometry of the streak line changes considerably com- 
pared to the co-rotating case as is illustrated in Fig. The 
stream line of the individual particle ejected from the stellar 
surface does not change and is ruled by the y6-law accelera- 
tion with the ambient wind in radial and by the conservation 
of angular momentum in azimuthal direction. However, since 
the spot from where the particles are ejected either rotates 
fixed relative to the stellar surface (vi-ot = Vs = -1-265 km s"') 
(Model 1) or moves in the opposite direction of the rotation 
(vi-ot = -H265kms"' and = -80.5 km s"^) (Model 2) the re- 
sulting streak line, i.e., the snapshot of all ejected particles for 
a given time, appears substantially different (cf. Fig. ll8t . 

Figure|5U|shows the correspondingly computed Rotbalmer 
model (Model 2). Despite being far from perfect, this model 
with a counter-rotating pulsation beat pattern and two streak 
lines emerging from the maxima of beat amplitudes does show 
some resemblance with the observed '8' or 'X' -shaped pattern 
in the dynamical spectra in Fig. |3 Particularly remarkable is 
the evolution of the pseudo-emission features from wind ve- 
locities beyond + v sin / to zero- velocity in nights 1 to 4 and the 
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Fig. 21. Temperature distribution on the stellar surface due to 
NRP beat patterns (Model 1) as seen from the equatorial plane 
at f = 48 hrs. The minimum and maximum temperatures in the 
visible regions are T = 21 500 K (black) and T = 25 000 K 
(white). 



phase of no variability in night 5, which are mostly responsi- 
ble for the resemblance with the observed 'crossing' features. 
Also the lack of features accelerating with the wind matches 
the observations better than the first model. Most important, 
also the time scales based on the beat period and the presence 
of two maxima in the beat amphtudes on the stellar equator 
match the observations surprisingly well. Clearly, the observed 
pseudo-emission and absorption features are more 'discrete' in 
velocity and time than the model prediction. But considering 
the large simplifications and the complete lack of physics in 
the transition zone from the pulsating photosphere to the base 
of the wind in our model, the result is encouraging for further 
detailed modeling. 

The important finding from our simplistic model therefore 
is, that the constructive interference of multi-periodic NRPs 
could play the crucial rule of providing the surface structure 
moving with the right time scales at the base of the wind to 
match (i) the observed base-of-the-wind variability (Ha) and 
(ii) the outer-wind structures rotating through the observer's 
line of sight as observed in the UV. 

7. Discussion 

The new high-resolution optical times series presented in this 
work allowed us to further investigate the line-profile variabil- 
ity of the B-supergiant HD 64760 in the photosphere and at the 
base of the wind. This new Feros data set provides detailed 
information to extend the Heros data set (KPS 2002). 

The high quality and the good sampling of the photospheric 
and wind variability over more than 2 rotational cycles of the 



star allowed us to reveal the presence of multi-periodic line- 
profile variations (LPVs) in the stellar photosphere. The three 
closely spaced periods could be identified as prograde high- 
order non-radial pulsation (NRP) modes with I and \m\ values 
in the range 6-8. The first detection o f high-order N RPs at a pe- 
riod of 0. 1 days had been reported by Baade ( 1984) - however, 
his data set was not adequate for a possible detection of closely 
spaced multi-periodic NRPs. Nevertheless, even the quality 
of our data set did not allow a very precise determination of 
the three pulsation modes from the small-amplitude (< 0.5%) 
LPVs. Therefore, we tested the two fundamental cases I - -m 
and I + -m via detailed modeling of the multi-periodic LPVs. 
Both cases allow us to represent the observations equally well. 
However, we favor the case / + -m to account for the fact that 
the observed 3% splitting into the three periods {P\ - 4.810hr, 
f 2 = 4.672 hr, P3 = 4.967 hr) is a true physical and asteroseis- 
mologically significant period splitting in the star's co-rotating 
frame, which implies that the m-values of the three modes have 
to be equal and different from /-values (our Model 2, cf. Tab.|3}. 
In the case of rotational splitting of pulsation modes, the ex- 
pected magnitude of the effect is such that the beating of ro- 
tationally split frequencies is of the order of the rotational fre- 
quency. Considering the measured split frequencies, it is found 
that their beat frequency is roughly consistent with the esti- 
mated rotational frequency of HD 64760. However, additional 
complications to this simple estimate have to be considered for 
an asteroseismological application of measured rotational split- 
ting ( Kawaler et al 1999) . It should be kept in mind that - if our 

preference for the / ^ -m case is correct - we are facing highly 
complex NRP surface patterns in HD 64760. For this reason we 
were restricted to use the geometrically simpler / - -m model 
for our subsequent modeling of stellar wind structure rooted in 
photospheric NRPs. 

The shortness of the measured periods and the absence of 
significant horizontal velocity fields clearly identifies the pulsa- 
tion modes as p-modes. The pulsation parameters with periods 
of a few hours as found for HD 64760 compare well wi th the 
findings for other O-type stars (cf. e.g. iHenrichsl (Il999l) for a 
compilation). It is particularly interesting that the first NRPs 
in early-type supe r giants were discovered in yArae (Blib) 
by iBaade & Ferletl ll 19841) with y Arae being almost a twin to 
HD 64760 with respect to the stellar parameters, including the 
very high value of v sin /. In particular the short-term variabil- 
ity in the photospheric lines of yArae with a derived period 
of 0.17 days and I - 10 is very reminiscent of the variability 
seen in this work in HD 64760. However, it is important to note 
that the variability in the UV wind-lines of yArae is quite dis- 
tinct f rom the periodic behavior of HD 64760, cf. IPrinia etall 

iFullerton et alJ (1 19971) had proposed as one possible source 
of regular surface structure in HD 64760 NRP patterns, which 
are locked to the 4.8-day rotation of the star, i.e., low-order 
g-modes with observed pulsation periods of 2.4 and 1.2 days 
and I - -m-values of 2 and 4, respectively. However, this sce- 
nario would imply an infinitely long pulsation period in the co- 
rotating frame (cf. Eqn.O, which appears not to be very phys- 
ical. Further, we have no evidence from models of pulsational 
instabilities for very long pulsation periods, for which the ro- 
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tation could become the dominant time scale in the observer's 
frame. 

KPS 2002 attributed the observed 2.4-day period in the Ha 
variability to modulations of the inner flanks of the emission 
humps at photospheric velocities to complex width variations 
of the underlying photospheric Ha profile and tried to connect 
these variations directly to the effect of NRP-induced LPVs of 
the photospheric Ha profile. The fact that the strongest vari- 
ability in the line profile is seen close to velocities of +v sin ; 
was taken as hint for the presence of strong horizontal velocity 
fields, i.e., the presence of low-order g-mode pulsations. The 
marginal detection of NRP-like LPVs in the He 1/14026 line fur- 
ther supported this hypothesis. With the new 2003 data set and 
the modeling results as presented in this work, this picture has 
to be refined. As shown in Sect. |5] NRPs are not expected to 
appear directly as measurable LPVs of the photospheric Ha 
profile due to smearing of the variations over the line profile 
as caused by the Stark-broadening of the intrinsic line profiles. 
Further, the new data set could not confirm the existence of 
NRP-induced LPVs on time scales of the wind modulation pe- 
riod but revealed the presence of multi-periodic NRPs, which 
are identified as p-mode pulsations with small horizontal ve- 
locity fields and periods of hours. Therefore, the observed Ha 
wind variations on time scales of days can only be attributed to 
secondary effects possibly linked to the presence of NRPs as 
discussed below. The interpretation of the observed Ha wind 
variations is further complicated by the presence of the dif- 
ferentially rotating wind: our simulations with the Rotbalmer 
code indicate that a simple dkect mapping of velocities in the 
line profile to a location in the stellar wind is not possible in the 
presence of rotationally twisted resonance zones. Indeed, any 
large-scale wind variability (like rotating streak lines) is pri- 
marily mapped into LPVs close to ±v sin The power distribu- 
tion over the line profile as e.g. derived from the models shown 
in Figs. ^91 and |20| display four power peaks blue- and red- 
wards of +vsin/ which nicely reproduces the observed power 
distribution of the observed Ha variability (cf. KPS 2002, their 
Fig. 3, blue and red peaks 'A' and 'D'). 

The most striking consequence from our finding of multi- 
mode NRPs with periods of a few hours is the presence of beat 
periods of several days. These beat periods can be directly seen 
in the observed and modeled photospheric LPVs as retrograde 
patterns traveling slowly (days) over the line profile superim- 
posed on the fast (hours) prograde patterns of the individual 
pulsation modes (see e.g. Fig.|8}. Simultaneously, we find the 
variability of the mostly wind-sensitive Ha line to be of the 
same time scales. Even some characteristics of the Ha line- 
profile variability like phases of low wind variability seem to 
coincide with phases of small beat amplitudes in the photo- 
spheric NRPs (cf. Fig.EJ. The He ii6678 line (cf. Fig.^, which 
shows on top of the photospheric NRP patterns some wind 
emission features, demonstrates directly the connection of the 
retrograde traveling NRP beat pattern with the wind emission, 
which apparently follows the beat pattern over the line profile 
(notably over the days 6-9). On the other hand, the NRP models 
of Ha have shown us that the photospheric NRPs themselves 
do not contribute to the LPVs of Ha (cf. Fig.fTTl. 



This is obviously an exciting finding since it provides us for 
first time with direct observational evidence for a connection 
of photospheric multi-periodic NRPs with spatially structured 
winds. The constructive or destructive interference of multi- 
ple photospheric pulsation modes with time scales of a few 
hours (the intrinsic stellar pulsation time scale of HD 64760 is 
< 1 day if a pulsation constant of Q - 0.04 days ( ILovv et all 
1984) and the stellar parameters in Tab. |2] are assumed) al- 
lows us to connect the fast photospheric variations at the base 
of the wind with the long wind modulation periods of days. 
Admittedly, the derived beat period of 6.8 days, which is ob- 
served in photosphere and base of the wind does not directly 
match the 1.2 and 2.4-day periods, which were derived for 
the periodic wind variability, but is more consistent with the 
longer 5.5 or 1 1-day repetition time scales as seen for the slow 
discrete absorption components (DACs) in the Iue Mega cam- 
paign data sets. In particular the simultaneous presence of the 
2.4-day period on one side and the wind emission features fol- 
lowing the slower beat pattern on the other side and the lack of 
an obvious link between the two periodicities inhibits a sim- 
ple unified picture of the photosphere - wind connection in 
HD 64760. Therefore, it appears particularly important to bet- 
ter understand the actual physical processes, which can connect 
photosphere structure to stellar wind structure. Without such an 
understanding it will not be possible to assess the significance 
of any observed and derived photosphere- and wind-variability 
periods. 

Fi gureOTIillustrates the spatial temperature structures in the 
photosphere created by the beating of the dominant two peri- 
ods in our NRP Model 1. The model predicts two diametrically 
opposite regions of maximum beat amplitudes (constructive in- 
terference) with temperatures ranging from T - 2 \ 500 K to 
T - 25 000 K at the stellar equator In the regions of destruc- 
tive interference 90 degrees away in azimuth, the temperature 
amplitudes are reduced and range from 22 300 K to 23 800 K. 
In our simplified scenario these NRP-induced spatial structures 
locally change the lower boundary conditions of the stellar 
wind, which in turn causes localized structural changes in the 
wind. (Note, that we used a simple 2 bright - 2 dark spot ge- 
ometry to represent the regions of constructive and destructive 
interference of NRP modes). 

The creation of wind structure by variation of lower bound- 
ary conditions have been investig a ted w ith hydrodynamical 
simulations bv lCranmer & Owockil ( 1 19961) who confirmed the 
formation of spiral structures (co-rotating interaction regions, 
CIR) originating from the collision of fast and slow streams 
rooted in the stellar surface. These wind structures modulate 
the observed stellar wind profiles as they are dragged through 
the line of sight to a distant observer by the underlying rotation 
of the stellar wind. 

In contrast to earlier models for the localized regions on 
the stellar surface of HD 64760, we investigated the effects of 
the observed and predicted retrograde motion of the regions 
of maximum beat amplitudes (equal to the regions causing the 
variation of the lower wind boundary condition) on the geome- 
try of the spiral structures and the resulting Ha LPVs. Despite 
the fact that our Ha simplistic model does not include any hy- 
drodynamics but is mostly geometrical, some key characteris- 
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tics of the observed Ha variations could be reproduced to some 
extent, which was not possible with models with the surface 
structures locked to the stellar rotation: the relevant time scales, 
and the velocity evolution of the pseudo-emission and absorp- 
tion features. 

The crucial question to be answered is whether the NRP- 
induced variations of the physical parameters like local temper- 
ature, gravity, and velocity fields are actually sufliciently large 
to introduce significant variations at the base of the wind to 
build up wind structures. ICranmer & Owocki(.1996) define the 
spot brightness amplitude A (with their Tq corresponding to our 



- equator 



) as: 



pot 



equator// equator 



Using the 'spot' temperatures from our NRP Model 1, we ob- 
tain amphtudes of 

A+ = ((25 000K)'*-(23 300K)4)/(23 300K)4 
= -1-0.325 

A- = ((21500K)'*-(23 300K)'*)/(23 300K/ 
= -0.275 

which are comparable to the A - +0.5 standard spot model 
used bv lCranmer & Owockil ill 9961) to model the Iue wind pro- 
file modulations of HD 64760 with CIRs. The azimuthal exten- 
sion of the standard spot model is with = 20 deg comparable 
to the azimuthal extension of the bright and dark sectors of our 
NRP Model 1. For the same model but with the temperature 
effects reduced to 50% to account for non-adiabatic effects (cf. 
Sect. 15. 4> we obtain spot brightne ss amplitudes of = 0.164 
and A" = -0.146 - according to 'Cranm er & Owocki ( Il996t) 
still sufficiently large to yield visible wind structure. 

While NRPs are not obviously very efficient in providing 
localized bright spots on the stellar surface for "Radiational 
Driven Orbital Mass Ejection" (RDOME) or "Pulsational 
Driven Orbital Mass Ejecti on" (PDOME) in the context of the 
creation of Be-star disks fowocki' 2004), they might be the 
source of the structured winds of hot supergiants like HD 64760 
as we have tried to demonstrate in this work. 

However, a hydrodynamical modeling, which takes fully 
into account the geometrical and physical boundary conditions 
as created by multi-periodic NRPs in the photosphere (as e.g. 
shown for the simple / - -m case in Fig. l21l i and which further 
includes the retrograde motion of the beat pattern with respect 
to the stellar rotation is still urgently needed to provide a con- 
sistent and time dependent picture of the stellar photosphere- 
to-wind structure of HD 64760. 
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